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SOLUTIONS OF THE LAMINAR COMPRESSIBLE BOUNDARY -LAYER 


EQUATIONS WITH TRANSPIRAT-ON WHICH ARE i 
APPLICABLE TO THE STAGNATION REGIONS 
OF AXISYMMETRIC BLUNT BODIES { 


By John T. Howe and William A. Mersman 


SUMMARY 


t 

‘Solutions of the laminar compressible boundary-layer equations én 
the ueu-end-Densughe=Sewm are obtained for Euler number 1/3 and Prandtl 
number at the wall equal to 0.7. Fifteen solutions are presented for 
various wall-temperature levels and transpiration rates. The solutions 
are applicable to the stagnation region of axisymmetric bodies by means 
of the Mangler transformation. UN 


INTRODUCTION . 


A detailed knowledge of velocity and temperature profiles in the 
stagnation region of blunt axisymmetric bodies having transpiration is 
of considerable interest. These profiles can be obtained from two- 
dimensional wedge flow solutions, in which the Ruler number is 1/3, by 
means of the Mangler transformation. Existing literature presents two- 
dimensional solutions for Eu.er numbers smaller and larger than this, 
but such that interpolation to obtain Euler number 1/3 solutions is con- 
sidered to be inaccurate. In order to cbtain accurate profiles for this 
Euler number, the equations of the two-dimensional laminar compressible 
boundary layer in the Brown, Donoughe, and Livingood form were integrated 
numerically. Solutions were obtained for five wall-temperature levels 
(Te/Ty = 0.25, 0.5, 1, 2, and 4), each having three transpiration 
conditions (fy = 0, -0.5, and -1.0). 


The purpose of this report is to fill a gap in the existing 
literature. Therefore, only brief discussions of the theory and method 
are presented, and the complete numerical results are listed. Transfor- 
mations to the desired axisymmetric coordinates are given so that local 
shear stress, heat transfer (as well as their corresponding dimensionless 
parameters), blowing rates, velocity profiles, and temperature profiles 
near the stagnation region of a blunt-nosed axisymmetric body can be 
obtained from the tabulated solutions. 
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Nu 


Re 


SYMBOLS 


constant used in equation (18) 
an arbitrary constant 


ety 


2 
Pelle 


local skin-friction coefficie.t, 


coefficient of specific heet at constant pressure 


dp/dx . 
Euler number, -x ——»s 
Pete 
dimensionless stream function aiid 
fe) c ; 
N oe 
heat-transfer coefficient, aw 
w- Te 
thermal conductivity ‘ 
an arbitrary reference length 
Nusselt number, hx 
ky 
Te-T 
temperature parameter, ane 9 oo 
Ty Ty 


ce 
Prandtl number, = 


heat-transfer rate, -k or 


Pyrex 
by 


Reynolds number, 


absolute temperature 

velocity parallel to the wall 

velocity normal to the wall 

distance along wall surface from stagnation point 


distance normal to wall surface : 


exponent of temperature in specific-heat relationship 


a 
€ exponent, of temperature in thermal conductivity relationship 
n independent variable in transformed boundary-layer equations, 
Wowe/ Myx 
T = Ty 
6 temperature variable, a 
e~ +w 
ve coefficient of viscosity 
p density 
T shear stress, yp ou 
oy 
wW exponent of temperature in viscosity relationship 
v stream function defined by equation (9) 
Subscripts 
e local conditions at the outer edge of the boundary layer 
w conditions at the wall 
Superscripts 
- axisymmetric conditions 


',",'" derivatives with respect to n 


THEORY 


The Differential Equations 


The two-dimensional partial differential equations describing the 
physical behavior of the leminar compressible boundary layer are: 


2 (pu) + = (pv) = 0 (1) 
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These equations comprise statements of conservation of macs, the momentum 


theorem, and conservation of energy, respectively. 
The boundary conditions are 
at y=0, 


u=0,vevif(x),T=T, (4) 


T+Te ,u Ue (5) 


The partial differential equations (1), (2), and (3) can be trans- 
formed to ordinary differential equations having only one independent 
variable, n, by use of the following relationships introduced by Blasius 
(ref. 1, p. 135), von Karman and Tsien (ref. 2), and Pohlhausen (ref. 1, 
p. 624), respectively. 


0,.u 
n=y TC (6) 
oy = ffiagpy (7) 
and . 
T-T 
6 xz = 
TN sa 


The stream function, y, in equation (7) is defined so that the continuity 
equation (1) is automatically satisfied; that is, 


a a 
and (9) 
2(0,,¥) 
eae. Ta 
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After certain assumptions and restrictions which vill be presented 
in the next section have been imposed, the two-dimensional partial dif- 
ferential equations, transformed to ordinary differential equations in 
the manner of Brow, Donoughe, and Livingood (refs. 3, 4, and 5), become 


plore 


Te 
L(g = (2) or + £°0' + (w+2)0'f" + (ee - weer 
| (10) 
and 
[Cero] on 


The complete derivation of the equations can be found in reference 3. 


The boundary conditions for equations (10) and (11) are as follows: 


at yn = O, 
f' =0 (12) 
f = f,, (required to be constant for similarity) (13) 
@=0 (14) 

at nwo 
@71 (15) 

T. 

ot f* - r (16) 


Boundary condition (12) is easily derived by means of the first of 
boundary conditions (4), the relationship being shown by equation (23). 
Boundary condition (13) is similarly derived by means of the second of 
boundary conditions (4), the relationship being shown dy equation (25). 
It is seen from equation (25) that f, is a measure of the transpiration 
rate. Negative values of fy correspond to bl or transpiration. 
Boundary conditions (14) and (15) are obvious, and (16) is derived by 
means of the last of boundary conditions (5), the relationships being 
given by equations (23) and (15). 
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Assumptions and Restrictions 


The derivation of equations (10) and (11) included certain simplifying 
assumptions (ref. 3). These will be listed and then diecussed in turn. 


1. The local Mach number at the outer edge of the boundary layer is 
small. 


2. The Euler number is constant. 
3. The wall temperature is constant. 


4, Local fluid properties within the boundary layer are proportional 
to some power of the absolute temperature. 


The first assumption makes it possible to omit the dissipative terms 
(the last two terms in eq. (3)) in the energy equation (see p. 13, ref. 6) 
because these terms are proportional to Mach number squared. For this 
same reason, T. and p, are considered to be constant. 


The second assumption permits integration of the Euler equation 


due 
- ea =z Pelle "ax (17) 


resulting in 
Ue = axEu (18) 


Equation (18) and therefore the two-dimensional solutions reported here 
are of the wedge flow type (ref. 7, p. 128). 


A combination of the first and third assumptions allows T,/Ty to 
be treated as a constant in equations (10) and (11). 


Although assumption number 4 states a proportional relationship of 
properties to some power of the ebsolute temperature, the constants of 
proportionality never have to be specified. The exponents, however, must 
be specified in order to integrate the differential equations. The power 
law relationships are simply 


rot @) (19) 
=) . 


Ses oY (2°) 

Cpw T. — 
-1 

a a 


The description of the two-dimensional theory is terminated at this 
point. A detailed discussion can be found in reference 3. Attention is 
now directed to a brief discussion of the method of solution, which differs 


somewhat from that of references 3 through 5. 
Method of Solution and Specified Numerical Data 


Solutions to equations (10) and (11) were obtained by numerical 
integration using the Adams-Moulton technique (ref. 8, p. 200, eqs. 6.6.2). . 
This is a variable interval method which is especially well suited to this 8 
problem. The numerical solutions were computed on the IBM 704 electronic 
data processing machine. 


The boundary conditions (12), (13), and (14) are not sufficient for 
starting the numerical integration of equations (10) and (11). Values of 
6' and f"' are also needed at 1 = 0 (i.e., Oy’ and fy"). These have to ‘ 
be estimated for the first attempt to integrate equations (10) and (11). . 
As n~o, the integration should yield the outer edge boundary conditions 
(15) and (16). Failure to produce and maintain these values indicates 
that the estimates for 6,' and f,'’ were in error and new estimates must 


be made. 
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In order to facilitate selection of the correct 6,' and f,,"', three 
initial trial integrations of a given example were made as follows: the 
first trial used estimates of 6,' and f,'’. In the second integration, 
only 6,' was changed from the first estimate, and in the third integra- 
tion only f,'' was changed from the first estimate. The results of these 
three trial runs yielded values for the partial derivatives of 6 and f' 
at n =o with respect to f,'' and 6,'. Linear interpolation then 
yielded revised estimates of the starting values, f,"' and Gy’. The 
process was then repeated until the correct initial values were obtained 
thus achieving the proper boundary conditions at 7s. 


Numerical data used for integrating equations (10) and (11) are 
Pry = 0.70 
Ws 0.70 


€ = 0.85 
a2 0.19 
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These numbers are applicable to air in the temperature range 1060° R 
to 2860° R according to reference 3. All examples were computed for 
Euler number 1/3 for reasons which will be presented in the section 
"Physical Relationships for Axisymmetric Stagnation Regions.” Examples 
were solved for five wall-temperature levels (Te/Ty equal to 0.25, 0.5, 
1, 2, and 4), each at three different transpiration conditions (f, = 0, 
~0.5, and -1.0). 


The numerical results of the integrations appear in tebles I through 
through XV. Velocity, temperature, and mass-flow profiles appear in fig- 
ures 1 through 9. Heat-transfer and skin-friction parameter curves are 
shown in figures 10 and ll. 


Physical Relationships for Two-Dimensional Wedge Flows 


The following relationships are useful in obtaining values for the 
physical quantities of interest from the tabulated results. Note that 
the left-hand sides of both equations (23) and (24) are tabulated in 


tables I through XV. 
=e ap + (we - 1) | (ref. 4) (23) 


pu » Te 
Dele f T, (ref. 4) (24) 


(ov), = -(B) (ref. 3) (25) 
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26) 
tv" wet y foes . ages (27) 
ce rie = 2f," (ref. 4) (28) 


ee (2) = -by(Te - ny)oy' | oe (29) 


oe = 6y' (ref. 4) (30) 


It can be seen from equation (25) that because f, is specified 


constant, v, is required to vary as x(Bu- 1)/ a 


These two-dimensional relationships (eqs. (23) through (30)) applied 
to the tabulated solutions for Eu = 1/3 (tables I through XV) correspond 
to the two-dimensional flow over a wedge. 


Physical Relationships for Axisymmetric Stagnation Regions 


Relationships resulting from the Mangler transformation will now be 
presented (ref. 7, pp. 129 and 168) which make shese wedge flow solutions 
applicable to the stagnation regions of bluni: axisymmetric bodies. This 
transformation is a stretching of the x and y coordinates so that the 
pressure, enthalpy, density, viscosity, and velocity parallel to the sur- 
face at the point x,y in the wedge flow are the same as those at X,y 
in the axisymmetric flow. In the following equations, quantities which 
apply only to the axisymmetric flow are designated with a bar. Those 
which apply to the wedge flow (as well as those mentioned above which 
apply to either flow) are written without a bar. The relationships derived 
by the Mangler transformation are as follows: 


aoe 
‘ate 


x = (ax) (31) 

¥-3y (32) 

oe a (33) 

Ki = 3m (34) 

WAZ) = H aylx) = - 43 wylte- Tey’ [= (35) 


Je} 


ile" 
36) 


+(%,F) -¢ +(x,y) = J3 uve fe [-+G- Ne] (- erors[(G- e+ 
( 
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Te a Mie - V5 op TY Vie = 3 ty" (37) 
Be = 3 Fe = V5 0 (38) 

(ov)y = F (ov)y = WB (B22), [a (39) 
wie = - = (40) 


Because the local Mach number exterior to the boundary layer is 
assumed to be small, the local exterior flow is considered to be almost 
incompressible. For the incompressible inviscid flow in the axisymmetric 
stagnation region, the velocity exterior to the boundary layer is a 
linear function of ¥ (ref. 7, p. 74, eq. (5.35)); that is, the axisym- 
metric Euler number (Bu) is unity in equations (34) and (39). This 
linear velocity for the axisymmetric stagnation region has also been 
observed experimentally over a wide range of free-stream Mach numbers. 
Equation (39) shows that the axisymmetric transpiration velocity, Vy, 
is constant because ue is proportional to ¥&. 


The two-dimensional Euler number (Eu) is 1/3 corresponding to Ru 
of unity according to equstion (34). All of the numerical solutions of 
equations (10) and (11) presented in this report are for Bu of 1/3 and 
can be used for the axisymmetric stagnation region. Applicecion of 
equations (31) through (40) to the solutions in tables I through XV will 
yield the desired axisymmetric blunt body stagnation region values. 
Profiles of velocity, mass flow, and temperature for the axisymmetric 
stagnation region are shown in figures 1 through 9 by using the appropri- 
ate ordinate. Heat-transfer and skin-friction parameter curves are 
presented in figures 10 and 11. 


DISCUSSION 


Some features of the solution and mumerical results varrant mention. 


The cases Te/Ty = 1 are somewhat special. The physical definition 
of @ (eq. (8)) reveals that either 6 is infinite or indeterminate for 
these exanples. However, the solutions for @ turn out to be finite, 
as is shown in figures 7, 8, and 9, and listed in tables IV, V, end VI. 
Therefore @, as defined physically, is actually indeterminate, that is, 
the temperature, T, is constant end is equal to T, (and T,) across the 
vouncary layer. However, @ is defined mathematically by the differential 


| 


eee 


uu 


equation (11) and the boundary conditions (14) and (15), but in this case 
it hay no physical meaning. For this degenerate situation, the sccond 
set of brackets in equation (10) is dropped, P is unity, and the equa- 
tion reduces to the familiar Falkner and Skan equation (ref. 1, p. 5.40) 
for constant property wedge flow conditions. 


Examination of figures 1 through 3 and tables I through VI shows a 
velocity overshoot condition for hot wall cases; that is, the velocity 
within the boundary layer exceeds that at the outer edge of the boundary 
layer. This results from the greater acceleration of the less dense hot 
gas in the boundary layer by a given (favorable) pressure gradient. It 
is seen that the amount of overshoot increases with increasing wall 
temperature or blowing rate. References 4 and 5 show that the amount of 
overshoot increases with Euler number. The two-dimensional flow at 
Bu «= 1/3 ses through the same values of u/u, as the axisymmetric 
flow at zl. Then it can be said that for a two-dimensional flow and 
an axisymmetric flow, both having Euler number unity, the magnitude of 
the velocity overshoot is smaller in the axisymmetric flow than it is in 
the two-dimensional flow. 


Fs os a 1, 5 Rew a: ug aie « es . i “a 
= As tan ee OE ee RY er 
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It is interesting to note in figures 4 through 6 and tables X 
through XV that an overshoot in mass flow occurs. Because of the increased 
density near a cold wall, this mass-flow overshoot is increased as the 
wall-temperature level is decreased, in contrast with the behavior of the 
velocity overshoot cond*tion. However, the amount of mass-flow overshoot 
increases with blowing rate as was the case with the velocity overshoot 


condition. ? 
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Figure 10 shows that for a given wall-temperature level, transpiration 
reduces the heat transferred to the wall to a remarkable extent. In par- 
ticular, for Te/Ty equal to 4, the heat-transfer parameter corresponding 
to transpiration parameter 0.667 is only about one fifth of its value for 
no transpiration. It is seen that the effectiveness of transpiration in 
reducing the heat-transfer parameter increases with the transpiretion rate 


and with 1T,/T,. 


As would be expected, figure 11 shows a similar effect of transpira- 
tion on skin friction. Quantitatively, for Te/Ty equal to 4, the skin- 
friction parameter corresponding to transpiration parameter 0.667 is only 
about one third of its value for no transpiration. 


SUMMARY OF RESULTS 


The Brown, Donoughe, and’ Livingood forms of the laminar compressible 
boundary-layer equations have been integrated numerically on the IBM 704 
electronic data processing mechine. Fifteen solutions were obtained for 
Euler number 1/3, wall Prendt] mmber 0.7, and for five ratios of local 
free-stream temperature to wall temperature each at three different 


12 


” 


transpiration conditions. These are T,/T,, values of 0.25, 0.5, 1, 2, 
and 4, and transpiration rates corresponding to fy of 0, -0.5, and -1.0. 
The solutions are presented in the forms of tables and graphs. 


Relationships are presented by which the tabulated solutions can be 
applied to stagnation regions of axisymmetric blunt bodies. 


The results indicate velocities in the boundary layer greater than 
the local exterior velocity for hot wal flow situations. The results 
also indicate local mass flows greater in the boundary layer than exterior 
to it for cold wall situations. These “overshoot” conditions are accentu- 
ated by increased transpiration for both the velocity profiles and the 
mass-flow profiles. 


A large reduction of skin friction and heat transfer is achieved by 
transpiration. The reduction is greater for situations in which the 
walls are colder than the local exterior flow. 


Ames Research Center 
National Aeronautics and Space Administration 
Moffett Field, Calif., Mar. 20, 1959 
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Figure 1.- Velocity profiles for f, = 0. 
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Figure 7.- Temperature profiles for fy = 0. 
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Figure 9.- Temperature profiles for fy = -1.0. 
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Fionre 10.- The influence of transpiration and wall temperature level on 
the heai-transfer parameter. 
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| Figure ll.- The influence of transpiration and wall temperature levei on 
i the ekin-friction parameter. 
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